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Summary. We compare all the available observational data on the redshift evo- 
lution of the total stellar mass and star formation rate density in the Universe 
with the mass and accretion rate density evolution of supermassive black holes, 
estimated from the hard X-ray selected luminosity function of quasars and active 
galactic nuclei (AGN). We find that on average black hole mass must have been 
higher at higher redshift for given spheroid stellar mass. Moreover, we find negative 
redshift evolution of the disk/irregulars to spheroid mass ratio. The total accretion 
efficiency is constrained to be between 0.06 and 0.12, depending on the exact value 
of the local SMBH mass density, and on the critical accretion rate below which 
radiatively inefficient accretion may take place. 



1 Introduction 

Observational evidence indicates that the mass of supermassive black holes 
(SMBH) is correlated with the luminosity (Marconi & Hunt, 2003 and refer- 
ences therein) and velocity dispersion (Trcmaine et al., 2002 and references 
therein) of the host spheroids, suggesting that the process that leads to the 
formation of galaxies must be intimately linked to the growth of the central 
SMBH. Studying low redshift AGN, Heckman et al. (2004) have shown that 
not only does star formation directly trace AGN activity, but also that the 
sites of SMBH growth must have shifted to smaller masses at lower redshift, 
thus mimicking the "cosmic downsizing" scenario first put forward to describe 
galaxy evolution by Cowie et al. (1996). Such a scenario has recently received 
many independent confirmations, both for the evolution of SMBH as traced 
directly by X-ray and radio luminosity functions (LF) of AGN (Marconi et al. 
2004; Merloni 2004; Hasinger et al. 2005), and for that of star forming galax- 
ies, thanks to large surveys such as SDSS, GDDS, COMBO-17, GOODS, etc. 
(see, e.g. Heavens et al. 2004; Juneau et al. 2005; Perez-Gonzalez et al. 2005; 
Feulner et al. 2005). 

Following Merloni, Rudnick and Di Matteo (2004; MRD04), here we dis- 
cuss a quantitative approach to the study of the posited link between star 
formation and SMBH growth, based on a detailed comparison of the redshift 
evolution of integral quantities, such as the total stellar mass, black hole mass 
and star formation rate densities. 
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Fig. 1. The upper left panel shows the evolution of the stellar mass density as 
a function of redshift, where the density is given as a ratio to the local value, 
p t .o = 5.6 x 1O 8 M0 Mpc -3 , the upper right panel shows the SFR density. Our best- 
fit model is also shown in each panel. Values of Ao = 0.3 and Pbh,o = 2.5 x 1O 5 M0 
Mpc -3 and m CI = are adopted here. Shaded areas represent 1-sigma confidence 
intervals of the model fits. The lower left panel shows a direct comparison between 
best-fit normalized mass density of spheroids (solid line, red shaded area) and black 
holes (dashed line, light red shaded area). Also shown is the evolution of the average 
Eddington scaled accretion rate rh(z) (dotted line, dark red shaded area). Finally, 
the lower right panel shows a direct comparison between the best-fit SFR and 
BHAR (rescaled by a factor ,4o(l + Ao)) densities. 



2 SMBH as tracers of galaxy evolution 

Under the standard assumption that black holes grow mainly by accretion, 
their cosmic evolution can be calculated from the luminosity function of AGN 
<f>(L\y \,z) = dN/dL^oi, where Lboi = eMc 2 is the bolometric luminosity 
produced by a SMBH accreting at a rate of M with a radiative efficiency 
e (Soltan 1982). Following the discussion in MRD04, we will assume that 
the absorption corrected 2-10 keV luminosity function of AGN, <fi(Lx, z), (La 
Franca et al. 2005) best describes the evolution of the entire accreting black 
holes population, yielding: 
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Lx is the X-ray luminosity in the rest-frame 2-10 keV band, and the bolomet- 
ric correction function Lboi(^x) is given by eq. (21) of Marconi et al. (2004). 
The exact shape of pnn(z) and ^bh^) then depends only on the local black 
holes mass density /?bh,o and on the (average) radiative efficiency e. This, in 
turn, is given by the product of the total accretion efficiency 77(a), itself a 
function of the inner boundary condition and thus of the black hole spin pa- 
rameter a, and a function / of the Eddington scaled dimensionless accretion 
rate m = iboi/^Edd- Below a critical rate, m cr , accretion does not proceed 
in the standard optically thick, geometrically thin fashion for which e=rj. 
Its radiative efficiency, instead, critically depends on the nature of the flow: 
if powerful outflows/jets are capable of removing the excess energy which is 
not radiated, as, for example, in the ADIOS scenario (Blandford & Begelman 
1999), then /=1, and black hole are always efficient radiators with respect to 
the accreted mass ("black holes are green!", Blandford 2005). On the other 
hand, if advection across the event horizon is the dominant process by which 
energy is disposed of (ADAF, Narayan & Yi 1995), we have: 



We use the total BHAR and mass densities, #bh(#) and pbh(-z) respectively, 
to estimate the redshift evolution of the average global Eddington scaled ac- 
cretion rate m(z) oc ^bh(^)/pbh(z) (lower left panel of Fig.^J and the corre- 
sponding radiative efficiency according to eq. ©. This allows us to identify 
the redshift at which a transition occurs in the global accretion mode of grow- 
ing SMBH. Depending on the assumed value of m cr , this transition redshift 
is 0, if black holes are always efficient accretors, i.e. m cr = 0, or z s=s 0.6 
if m cr rs 0.05. Given the overall evolution of the BHAR density, this also 
implies that radiatively inefficient accretion could contribute to only a small 
fraction of the total black hole mass density (see also Yu and Tremaine 2002; 
Merloni 2004; Hopkins et al. 2006 and references therein). 

2.1 The parallel evolution 

Our goal is to link the growth of SMBH from eq. to the growth of stel- 
lar mass in galaxies. Because local SMBH are observed to correlate with 
spheroids only, we introduce the parameter A(z), the ratio of the mass in disks 
and irregulars to that in spheroids at any redshift, so that the total stellar 
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Fig. 2. 1,2 and 3 sigma confidence contours for the average accretion efficiency 
r] and the index a describing the evolution of the stellar spheroid to black hole 
mass density ratio. In each panel, the leftmost and rightmost set of contours show 
the results obtained assuming pBH.o = 4.5 (Marconi et al. 2004), and 2.5 x 1O 5 M0 
Mpc -3 (Yu and Tremaine 2002), respectively. Also shown is the dependence of such 
constraints on the assumed value of the critical rate between radiatively efficient 
and inefficient accretion (ADAF-like), with the dashed contours in the right panel 
showing the m cr =0 reference point. 



mass density can be expressed as: p*(z) = p sp h{z) +/9disk+irr(z) = p sp h(z)[l + 
X(z)]. We then assume that X(z) evolves according to X(z) = A (l + z)~@ , 
where Ao is the value of the disk to spheroid ratio in the local universe. Also 
we assume that the mass density of spheroids and supermassive black holes 
evolve in parallel, modulo a factor (1 + z)~ a , obtaining a prediction for the 
observable stellar mass density evolution as traced by SMBH growth: 

p,(z) = Ao P Bn(e, z)(l + z)- a [l + A (l + z) - ' 3 ] (4) 

where Ao is the constant of proportionality in the Magorrian relation. By 
taking the derivative of (QJ, accounting for stellar mass loss, an expression is 
also found for the corresponding star formation rate (SFR) density evolution 
(see eq. (7) of MRD04). 

With these expressions we obtain statistically acceptable simultaneous 
fits to all available observational data points (see MRD04 for a complete list 
of references) of both p*(z) and SFR(z). For each choice of pbh,o, Ao, and 
of the critical accretion rate m cr , the fitting functions depend only on three 
parameters: a, f3 and the accretion efficiency r/. One example of such fits is 
shown in Fig.Q]for the specific case Pbh,o = 2.5 x 10 5 M© Mpc -3 , Ao = 0.3 and 
m a = 0. Because the drop in the AGN integrated luminosity density at low z 
is apparently faster than that in SFR density (see lower right panel of Fig. , 
the average black hole to spheroid mass ratio must evolve with lookback time 
(a > 0; see lower left panel of Fig. [JJ. This result is independent from the 
local black hole mass density, or from Ao , and is not strongly affected by the 
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choice of the value for the critical accretion rate. This is shown also in Fig. El 
where the constraints on the fit parameters are shown as confidence contours 
in the accretion efficiency-a plane for various choices of m cr . As SMBH grow 
most of their mass at high m, the global constraints on rj are not strongly 
affected by any reasonable choice of m C r- 

3 Conclusions 

We have made quantitative comparisons between the redshift evolution of the 
integrated stellar mass in galaxies and the mass density of SMBH. Although 
clearly correlated with the stellar mass density, SMBH accretion does not 
exactly track either the spheroid nor the total star assembly: irrespective 
of the exact mass budget in spheroids and disks + irregulars, the ratio of 
the total or spheroid stellar to black hole mass density was lower at higher 
redshift. Our results also suggest that the fraction of stars locked up into 
the non-spheroidal components of galaxies and in irregular galaxies should 
increase with increasing redshift ((3 < 0, see MRD04 for a discussion of this 
point). Our version of the Soltan argument yields a well defined constraint 
on the average radiative efficiency, and a corresponding one for the total 
accretion efficiency (i.e. on the mean mass- weighted SMBH spin), only weakly 
dependent on the uncertain physics of low m accretion flows. 
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